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The induction equation of kinematic magnetohydrodynamics is mathematically equivalent to 
a system of integral equations for the magnetic field in the bulk of the fluid and for the electric 
potential at its boundary. We summarize the recent developments concerning the numerical imple- 
mentation of this scheme and its applications to various forward and inverse problems in dynamo 
theory and applied MHD. 
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1 Introduction 

This paper is intended as a summary and update of our recent work on the theo- 
retical formulation, the numerical implementation, and the practical application of 
the integral equation approach (lEA) to kinematic magnetohydrodynamics (MHD). 
Actually, our activity in this field had started with the paper (Stefani et ai, 2000), 
written together with Karl- Heinz Radler, whom we would like to thank not only for 
his insistence on a rigorous mathematical formulation in this particular case, but 
also for many fruitful discussions over the course of time. 

The theory of hydromagnetic dynamos is mainly concerned with the self-excitation 
of cosmic magnetic fields, with particular focus on the fields of planets, stars, and 
galaxies (Krause and Radler 1980). As long as the self-excited magnetic field is 
weak and its influence on the velocity field is negligible we speak about the kine- 
matic dynamo regime. When the magnetic field has grown to higher amplitudes 
the field-generating velocity field is modified, and the dynamo enters its saturation 
regime, which will not be considered in the present paper, though. 

The familiar way to deal with kinematic dynamo action relies on the induction 
equation for the magnetic field B, 



where u is considered as a pre-given velocity field, /xq is the magnetic permeability 
constant, and a is electrical conductivity of the fluid. The time dependence of the 
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magnetic field B in Eq. ([T]), in terms of exponential growth or decay, is governed 
by the ratio of field changes due to velocity gradients (first term on the r.h.s.) 
to the field dissipation (second term on the r.h.s.). This ratio is determined by 
the magnetic Reynolds number = fJ^cxLU, where L and U are typical length 
and velocity scales of the flow, respectively. When the magnetic Reynolds number 
reaches a critical value, denoted by i?^, the field can grow exponentially in time. 

Equation ([1]) follows directly from pre-Maxwell's equations and Ohm's law in 
moving conductors. In order to make this equation solvable, boundary conditions 
of the magnetic field must be prescribed. In the case of vanishing excitations of the 
magnetic field from outside the considered finite region, the boundary condition of 
the magnetic field is given by B = 0(r~^) as r — oo. For dynamos in spherical 
domains, as they are typical for planets and stars, the problem of implementing the 
non-local boundary conditions for the magnetic field is easily solved by using de- 
coupled boundary conditions for each degree of the spherical harmonics. For other 
than spherically shaped dynamos, in particular for the recent laboratory dynamos 
working often in cylindrical domains, but also for dynamos in galaxies and plan- 
etesimals, the handling of the non-local boundary conditions is a notorious problem 
(Stefani et al. 2009). 

A simple, but rather rude way to circumvent this problem is to replace the com- 
plicated non-local boundary conditions by simplified local ones (so-called "vertical 
field conditions"). This method is sometimes used in the simulation of galactic 
dynamos (a recent example is Hubbard and Brandenburg 2010). 

For the simulation of the cylindrical Karlsruhe dynamo experiment, Radler et 
al. (1998, 2002) had used the alternative trick of embedding the actual electrically 
conducting dynamo domain into a sphere, and the region between this sphere and 
the surface of the dynamo was virtually filled by a medium of lower electrical con- 
ductivity. By reducing the value of this conductivity, it was possible to check the 
convergence of the numerical solution. 

Of course, both methods are connected with losses of accuracy, the latter method 
being certainly more accurate. In order to fully implement the non-local boundary 
condition, Maxwell's equations must be fulfilled in the exterior, too. This can be 
implemented in different ways. For the finite difference simulation of the Riga dy- 
namo, a Laplace equation was solved, for each time-step, by a pseudo-relaxation 
method, in the exteriour of the dynamo domain and the magnetic field solutions in 
the interiour and in the exterior were matched using interface conditions (Stefani et 
al. 1999, Kenjeres et al. 2006). A similar method, although based on the finite ele- 
ment method, was developed by Guermond et al. (2003, 2007) and was intensively 
used later for practical dynamo applications (Guermond et al, 2009, Giesecke et al. 
2010a, b, Nore et al. 2011). 

A quite elegant technique to circumvent the solution in the exteriour was pre- 
sented by Iskakov et al. (2004, 2005) who used a combination of a finite volume and 
a boundary element method. This latter method was recently used by Giesecke et al. 
(2008, 2010a, b) for the simulation of the French VKS-dynamo (although the main 
focus of that work laid on the particular effect of the high magnetic permeability 
impellers on the mode selection in this type of dynamo). 

An alternative to the solution of the induction equation is the integral equation 
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approach (lEA) for kinematic dynamos which basically relies on the self-consistent 
treatment of Biot-Savart's law. For steady dynamo action in infinite domains of 
homogeneous conductivity, the lEA becomes very simple and had already been 
employed by a few authors (Gailitis 1970, Gaihtis and Preibergs 1974, Dobler and 
Radler 1998). For dynamo action in finite domains, however, the simple Biot-Savart 
equation has to be supplemented by a boundary integral equation for the electric 
potential (Roberts 1967, Stefani et al. 2000, Xu et al. 2004a). If the magnetic field 
becomes time-dependent, yet another equation for the vector potential can be added 
to get a closed system of integral equations (Xu et al. 2004b). 

In the following section, we will present this formulation and its specifications in 
detail. Then we will consider its application to various dynamo problems, including 
those with different degrees of symmetry which enable dimensional reductions. At 
last, we will go over from forward problems to inverse problems for which the lEA 
represents an ideal mathematical starting point. 



2 Mathematical formulation 

Assume a fluid of electrical conductivity a, flowing with velocity u, to be confined in 
a finite region V with boundary S", the exterior of this region consisting of insulating 
material (or vacuum). By further assuming the velocity u field to be stationary, we 
may write the electric field, the magnetic field and the magnetic vector potential in 
the following form: 

E(r, t) = E(r) exp At, B(r, t) = B(r) exp At, A(r, t) = A(r) exp At, (2) 

where the real part of A is the growth rate, and its imaginary part is the frequency of 
the fields. In the general case, we allow external currents to be present that produce 
a magnetic field Bq. The remaining magnetic field, induced by the flow, is denoted 
by b. Then, the current density is governed by Ohm's law 

j = a(E + ux (Bo + b)) 

= a{-V(p - AA + u X (Bo + b)) (3) 

where the electric field E, the electric potential (j) and the magnetic vector potential 
A are used. Applying now Biot-Savart's law to the magnetic field. Green's sec- 
ond theorem to the solution of the Poisson equation for the electric potential, and 
Helmholtz's theorem for the integral representation of the vector potential in terms 
of the magnetic field, we arrive at the following integral equation system: 

b^r) ^ f (u(rOx(Bo(rO + b(rO))x(r-rO ^^^, 
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Afr) = - f (Bo(^0 + b(rO)x(r-rO ^^^, 
4717^ |r-r'|3 

/ Bo(sO + bM^^, 

where /x is the permeabihty of the fluid. 

Note that this integral equation system is by far not the only possible one. The 
double use of the magnetic field and its vector potential might even seem a bit 
awkward. There are indeed other possible schemes, one of them starting from the 
Helmholtz equation for the vector potential which leads, however, to a nonlinear 
eigenvalue problem in A (see Dobler and Radler 1998), while the above scheme is a 
linear eigenvalue in A which has many advantages when it comes to the numerical 
treatment. 

The general integral equation system (jUlH]) can now be further specified. In the 
case with Bq 7^ and small Rm, i-e. below the threshold of self-excitation, the 
system describes an induction problem in which the applied field Bq is only slightly 
deformed by the velocity. In the case of Bq = 0, the system represents an eigenvalue 
problem for the complex constant A, with only negative eigenvalues for small Rm 
(below the dynamo threshold), and one or a few eigenvalues with positive real parts 
for the case of large R^. (above the dynamo threshold). 

In the case of a steady problem, with neither any growth/decay nor any oscilla- 
tion of the field, this equation system reduces to 



^ I (u(rOx(Bo(rO + b(rO))x(r-rO ^^^, 

47r ]y |r - r'l3 



1 _ 1 /- (u(r')x(Bo(r') + b(r')))■(s-r') 



^^V^' 



2 A-K jy |s - r'P 

Again, we can distinguish the case with Bq 7^ and small Rm for which it repre- 
sent a stationary induction problem. This is of particular interest for many technical 
liquid metal problems (steel casting, crystal growth) characterized by a rather small 
Rm- In this case, the induced field b can be omitted under the integrals which leads 
to a linear relation between the source u and the result b of the induction process. 
This linear relationship represents a convenient starting point for the treatment of 
the inverse problem to infer u in the bulk of the fluid from values of b measured in 
the exteriour of the fluid. A problem of this sort will be discussed in the next but 
one section. In the case of Bq = 0, the steady system represents again an eigenvalue 
problem, but now for the critical Rm that leads to a marginal and non-oscillatory 
dynamo. 
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3 Numerical implementation and dimensional re- 
ductions 



In this section we will illustrate the general method by treating dynamo problems 
first in 3D and then in 2D and ID, the latter cases being characterized by different 
degrees of symmetries that can be exploited for a dimensional reduction of the 
problem. We will see that, in stark contrast to the corresponding procedure in the 
differential equation approach, such a dimensional reduction represents a formidable 
task in the IE A. 

3.1 3D - Matchbox dynamos 

Let us start with a full problem in 3D for which we will summarize the result of 
the papers (Xu et al. 2004a, 2004b) for a so-called "matchbox dynamo", i.e. a 
mean-field dynamo in a rectangular box. The original assumption of a full velocity 
field as the source of induction can be easily translated to the case with a helical 
turbulence parameter a. Assume that we apply a specific spatial discretization of 
all fields in Eqs. (jHE]), we formally obtain 

hi = fia[Pik{Bok + bk) - XRijAj - Qii(f)i], (9) 

Gml4>l = SmkiB^k + i>k) — XTmjAj, (10) 

Aj = WjkiBok + bk), (11) 

where Einstein's summation convention is assumed. We have used the notion Gmi = 
0.5 Smi + Umi- -Bofc and denote the degrees of freedom of the externally added 
magnetic field and the induced magnetic field, Aj the degrees of freedom of the 
vector potential in the volume V, (pi the degrees of freedom of the electric potential 
at the boundary surface. For the later numerical utilization with different induction 
sources it is useful that only the matrices Pik and Smk depend on the velocity (or 
on the corresponding mean-field induction sources), while Rij, Qu, T^j, Gmi and 
Wjk depend only on the geometry of the dynamo domain and on the details of the 
discretization. 

Inserting Eqs. (|TOi) and (ITTi) into Eq. ([9]), which eliminates Aj and (pi, we end 
up with one single matrix equation for the induced magnetic field components bf. 

bi = fi(T[Pik{BQk + bk) — \RijWjk{Bok + bk) — QiiGj^Smk{.BQk + bk) 

+\QiiGY^T^,Wjk{Bok + bk)] (12) 

which can be further transformed into 

[5ik - fiaEik - fia\Fik]bk = [fitrEik + fia\Fik]Bok, (13) 

where Eik = Pik - QuGj^Smk and Fik = -RijWjk + QnG^^TmjWjk- It should be 
noted, though, that the needed inversion of the boundary integral equation (|T0|) for 
the electric potential needs some specific care (Wondrak et al. 2009). This has to do 
with the singular character of this equation which, in turn, mirrors the ambiguity 
of the electric potential with respect to some additive constant. 



5 



For the computation of induction effects, the induced magnetic field can be 
obtained by solving the algebraic equation system f lT2|) . In order to compute kine- 
matic dynamos without external sources, Eq.f ll2p reduces to the following general- 
ized eigenvalue problem for A: 

[5ik - ficrEik]bk = XficrFikbk, (14) 

To test the working of the scheme we apply it now to dynamos in rectangular 
geometry, which we would like to coin "matchbox dynamos". As usual in mean- field 
dynamo theory (see Krause and Radler 1980), a parametrizes the induction effect 
of helical turbulence on the large scale magnetic field. Formally speaking, the term 
u X B in the electromotive force term in Eq. ([3]) has to be replaced by a ■ B. 

While the spherically symmetric dynamo has played a paradigmatic role in 
dynamo theory, its counterpart in rectangular geometry is certainly a highly artificial 
problem without any application in astrophysics. Nevertheless, it may illustrate the 
capabilities of the lEA to solve dynamo problems in geometries that are far away 
from spherical. 

In our specific example (for more cases see Xu et al. 2004b) we consider a 
homogeneous distribution, a(r) = C, in a "matchbox" with side lengths 2.0x1.6x1.2 
which gives approximately the same volume as a corresponding sphere of radius 1. 
The domain is divided into 10^ smaller rectangular boxes, with each face divided 
into 10^ rectangles. 



EF 1 EF 2 EF 3 




Figure 1: Magnetic eigenfields corresponding to the three dominant eigenvalues of 
the dynamo in a matchbox with sidelengthes ratio 2x1.6x1.2. 

For this setting, we find the first three critical values of C to be C'l = 4.728, 
6*2 = 4.898, and C3 = 4.934. The corresponding first three eigenfunctions are 
shown, for the free decay case C = and the slightly sub-critical case C = 4, in 
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figure 1. These eigenfunctions still resemble the corresponding eigenfunctions in the 
spherical case. However, the exact degeneracy of the eigenvalues in the latter case 
is lifted now. The lowest eigenvalue corresponds to an eigenfunction whose dipolar 
axis is perpendicular to the largest face. Nevertheless, the three eigenvalues are still 
close the well-known value 4.49 for the spherical case. 



3.2 2D - VKS-like dynamos 

The recent successful dynamo experiments (see Stefani et al. 2008 for a recent 
survey) in Riga (Gailitis et al. 2000), Karlsruhe (Sticglitz and MiiUer 2001), and 
Cadarache (Monchaux et al. 2009) were all carried out in cylindrical vessels filled 
with liquid sodium. For this reason it is worth to specify the integral equation ap- 
proach to this geometry. As long as the dynamo source (i.e. the velocity field or a 
corresponding mean-field quantity) is axisymmetric, the different azimuthal modes 
of the electromagnetic fields with the angular dependence according to exp(im</?) 
can be decoupled. Ultimately, this leads to a tremendous reduction of the numer- 
ical effort. The price we have to pay for this is the necessity to carefully deriving 
the dimensionally reduced version of the integral equation system. The necessary 
integrations over (p turnes out to be a formidable task (see Xu et al. 2008), quite 
in contrast to the respective procedure in the differential equation approach, where 
the expression d/dip can simply replaced by im. 

Lets assume now the electrically conducting fluid to be confined in a cylinder 
with radius R and total height 2H. Introducing the cylindrical coordinate system 
(p, (f, z), we have 

r = [p cos (f, psm(p,zf ,h = [bp, b^, bzf,u= [up, u^, u^Y . (15) 

The magnetic field b, the electric potential 0, and the vector potential A all are 
expanded into azimuthal modes: 




E 




exp(im</7). (16) 



As long as the velocity field is axisymmetric (i.e. it has only a component with 
m = 0), the fields [bj„,0^,A^]^ with different m — 0, ±1,±2, ••• decouple from 
each other. For the sake of convenience, we will always re-denote [b^, 0m; A^]-^ as 
[b, 0, A]^. The reduct ion to a problem exclusively in r and z is then achieved by 
carrying out all integrations over 0. This painstaking procedure results finally in 
the following system of matrix equations 
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A A ( + bp\ 

= W\ Bo^ + b^\, (19) 
A J \ Boz + bj 

where the matrix elements of P, Q, R, S, T, U, and W can be found in the paper 
by Xu et al. (2008). Combining Eqs.([I7}lIlD, we obtain 

(I - fiaE - fiaXF)h = fia{E + AF)Bo, (20) 

where 

E = P-Q-(^H-U)-^-S (21) 

F = Q-(^I + U)"^-T-W-R-W. (22) 

Again, induced magnetic fields can be obtained by solving the algebraic equation 
system fl2U]) . while for kinematic dynamo problems, the generalized eigenvalue prob- 
lem 

(I - fiaB) ■ b = fiaXF ■ b (23) 

has to be solved. 

In connection with the optimization of the VKS dynamo experiment, Marie, 
Normand and Daviaud (2006) had studied an analytical test flow of the same topo- 
logical type as the flow in the real experiment (flow topology s2+t2, what means 
two poloidal eddies with radial inflow in the equatorial plane, together with two 
counter- rotating toroidal eddies). The velocity field of this "MND flow" reads: 

Ur = -|r (l-r)2(l + 2r)cos(7rz) (24) 

= 4er(l - r) sin(7rz/2) (25) 
Uz = (1 — r)(l -|- r — 5r^) sin(7r2;). (26) 

For the parameter e, which determines the ratio of toroidal to poloidal flow, we have 
used the value e = 0.7259. For the case without any side layer, the structure of 
the magnetic eigenfield is illustrated in figure 2. In figure 2a the field lines of the 
equatorial dipole are seen, while the isosurfaces of the magnetic field energy (figure 
2b) show the typical banana-like structure. 

3.3 Reduction to ID - Spherically symmetric dynamo 

After having utilized, in the former subsection, the cylindrical symmetry of some 
flow field for reducing the full 3D lEA system to a 2D system, we go now one 
step further and invoke spherical symmetry to make a reduction to a ID problem. 
This is possible for a spherically symmetric a"^ dynamo. In contrast to the original, 
analytically solvable model with constant a (Krause and Radler 1980), we will allow 
the profile a to vary with the radial coordinate r. After summarizing the derivation 
of the two coupled radial integral equations as given in (Stefani et al. 2000, Xu et 
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(a) 



(b) 



Figure 2: Simulated magnetic field structure of the dominant eigenmode of the MND 
flow, (a) Magnetic field lines of the equatorial dipole field, (b) Iso-surface of the 
magnetic field energy. 



al. 2004b), we will treat numerically a non-trivial model with a radial dependence 

a{r) = Cr'^. 

The divergence- free magnetic field B is split, as usual in dynamo theory, into a 
poloidal and a toroidal part, denoted by Bp and By. Using the Coulomb gauge, 
V • A = 0, the same can be done for the vector potential, A = Ap + At- All these 
fields can be represented by the defining scalars S, T, 5"^, according to 



Bp = V X V X ( -r 



B. 



V X ( -r 

r 



(27) 



Ap = VxVx( — rV At = Vx (— r 
r / \ r 



(28) 



For our spherical problem, the defining scalars and the electric potential can be 
expanded in series of spherical harmonics Yim{0, 0). As an example we indicate 



S{r,9,<l>) - J]sUr)l'/U^» , 

l,m 



(29) 



and corresponding expressions can be written for T{r,9,(f)), S^{r,9,(l)), T^{r,9,(l)) 
and ip{r,9,(f)), with sim{r) being replaced by tim{r), s^(r), tf^{r), and ipim{r), re- 
spectively. 
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For the spherical harmonics Yj 



Im \ 



the definition 



Y, 



'2/ + 1 (/ -m)! 



Pim(cose)e" 



irruj) 



Att {l + m)\ 

is employed, which implies the following orthogonality relation for the Yim{ 



(30) 



2n 



sine d9 Yit^, {9, 4>)Yi 



ll'^mm' 



From Eqs. 



we obtain the components of B in the form 

l,m 

^ 1^ timjr) dYi^{e,(j)) ^ 1 dsi^{r) dYim{9, 0) 



Lm 



r sin^ 



tirair) dYi 



lm \ 



dO 
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r dr do 

1 ciszm(r) dYi^ 



r sin 6' 



(31) 

(32) 
(33) 
(34) 



and equivalent expressions for the components of A, in which sim{T) and tim{r) are 
replaced by s^(r) and t^{r), respectively. 

A necessary ingredient for the dimensional reduction is the expression for the 
inverse distance between two points r and r' in terms of spherical harmonics, 



oo I 



- + IrL 

1=0 m=-l > 



(35) 



where r> denotes the larger of the values r and r', and r< the smaller one. 

On this basis, the two coupled integral equations for the functions sim{r) and 
Umij) can be derived. The first equation for s/m(r) can easily be obtained by multi- 
plying the magnetic field with the unit vector in radial direction and then integrating 
over the angles 9 and 0. The derivation of the equation for tim{r) needs more work, 
including the treatment of the electric potential and the vector potential at the 
boundary. 

Here we go straight to the final result of this procedure, whose details were 
worked out by Xu et al (2004b): Sim{r) and tim{r): 



21 + 1 



air') t,„Ar') dr' + 



a{r')tiUr') dr' 



ji+i 



X / — —sim{r')dr' - X 



sim{r')dr' 



(36) 



and 

Umir) 



a{r)sim{r) 



21 



da{r' 



—sim{r')—dr' 



1 ./o d'^' ' ' 



10 



R 



21 + 1 dr 



da{r') r'+^ / + 1 r'+^ 



2/ + 1 i?2i+i 
A 



r tira{r )dr 



21 + 1 Jq r 



-^ti^{r')dr' 



21 + \ 



tim{r')dr' - -^^a{R)sim{R)] 



(37) 



The correctness of this system can be easily checked by differentiating equations 
(l36l) and (!37j) two times with respect to the radius which gives the two well-known 
differential equations for Simij) and Umir): 



d'^Slm _ l(J + 1) 

dr"^ r 



-Sir 



+ a{r)tir. 
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dr"^ ^2 



tlr. 



d_ / , . (jg;^ 

dr \ dr 



l{l + l) , , 
H 5 — OL[r)sir, 



(39) 



with the boundary conditions 



tlmiR) = R 



dsim{r) 



dr 



+ IsUR) = 0. 



(40) 



It should be noted that in this particular case the complicated dimensional re- 
duction (by expressing the fields in spherical harmonics and integrating over the 
spherical angles) can be circumvented by a simpler one. This starts directly from 
the differential equation system fl38H40l) for which the Green functions can be derived 
which then lead to the integral equation system (13611371) (Xu et a/.2004b). 

Now we turn to the numerical utilization of the radial integral equation system 
(13611371) . The discretization in radial direction was done with a grid number of 
128 which gives already a satisfactory accuracy. The eigenvalue solver gives us 
immediately all the eigenvalues. This is illustrated in figure 3, which shows for the 
case a(r) = Cr^ the real part (figure 3a) and the imaginary part (figure 3b) of 
the first 8 eigenvalues. The existence of two so-called exceptional points, at which 
two real eigenvalue branches coalesce and continue as a pair of complex conjugated 
eigenvalues, is clearly seen. 



4 Inverse problems for small 

Assuming a velocity field as given, and asking for the resulting magnetic field, rep- 
resents a typical forward problem. The question can, however, also be inverted: 
assume the magnetic field as being measured in the exteriour of the fluid, what is 
the velocity that produces this field? This is a typical inverse problem, and it has 
puzzled dynamo theorist since many years. 

In the general case, for arbitrary values of Rm-, it represents a highly nonlinear 
inverse problem. Some progress in its treatment has been made for restricted set- 
ups, for example by applying the so-called frozen-flux approximation for the Earth's 
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Figure 3: Eigenvalues of the spherically symmetric a dynamo with a radial depen- 
dence a{r) = Cr^. (a) Growth rates, (b) Frequencies. 

core which allows to determine (still with some appropriate regularization) solutions 
for the tangential velocity field components at the Core-Mantel boundary from the 
measured time dependence of the radial magnetic field components (Roberts and 
Scott 1965). A complementary sort of restricted models is concerned with the deter- 
mination of the radial dependence of a{r) of an assumed spherically symmetric a"^ 
model from spectral properties. Solving such a type of restricted inverse problems 
(by means of an Evolutionary Strategy) it was possible, e.g., to obtain such a pro- 
files that lead to oscillatory dynamo solutions (Stefani and Gerbeth 2003). Apart 
from those special solutions, inverse dynamo theory is still in its infancy. 

The situation becomes much better for the case of induction problems with small 
Rm- As mentioned above the induced field b can then be omitted under the integrals, 
and we obtain a linear relation between the (wanted) source u and the (measured) 
result b of the induction process. 

Although it is far from being a generic geo- or astrophysical problem, we would 
like to illustrate the typical solution procedure of such inverse problems by a model 
related to industrial steel casting. Figure 4a shows the set-up of a corresponding 
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model experiment in which a hquid metal (in our case GalnSn) is poured from a 
tundish through an submerged entry nozzle into the mould. Among other problems 
of industrial interest, like flow instabilities due to Argon entrainment into the metal, 
we deal here with the problem how a magnetic stirrer can influence the flow structure 
in the mould. For this purpose, we use a so-called one-port submerged entry nozzle 
with just a hole at the bottom of the nozzle. 

Based on the integral equation system ([3IH]), specified to small Rm, we have 
developed the so-called Contactless Inductive Flow Tomography (GIFT) for the 
reconstruction of velocity fields from externally measured magnetic fields (Stefani 
and Gerbeth 1999, 2000a, b, Stefani et al. 2004). While GIFT is in principle able to 
infer full 3D velocity fields by applying subsequently two different (e.g. orthogonal) 
external magnetic fields, for the present case of thin slab casting it can be reduced 
to the determination of the velocity components parallel to the wide faces of the 
mould (Wondrak et al. 2010). For this it is enough to apply only one magnetic 
field by a single coil (see figure 4a). The interaction of the flow with the applied 
field produces induced magnetic fields that we measure at a number of positions at 
the narrow faces of the mould in order to reconstruct from them the velocity field. 
The mathematics of this inversion relies in the minimization of the mean squared 
deviation of the measured magnetic fields from the fields resulting according to the 
integral equation system (I71|8]) from the velocity field. This minimization is done by 
solving the normal equations, whereby we use various auxiliary functionals which 
serve to ensure the divergence-free condition of the velocity, to enforce its two- 
dimensionality, and to minimize in parallel its mean quadratic value, weighted by 
some regularization parameter (Tikhonov regularization). Figure 4b shows a time 
sequence of reconstructed velocity fields resulting from this inversion, for a particular 
flow experiment with applied stirring. What is clearly seen from the four plots is 
a significant (left-right asymmetric) up- and down movement of the vortex centers 
(indicated by the dark patches). 

5 Conclusions 

In this paper we have surveyed the principles and some applications of the integral 
equation approach (lEA) to kinematic MHD. The lEA has turned out a viable 
scheme for the correct numerical treatment of dynamo problems in non-spherical 
domains, for examples in "matchboxes" and, more importantly, in cylinders (see 
also Avalos-Zuniga et al. 2007). It has also proved a good starting point for the 
treatment of inverse problems of MHD, with direct applications in a number of 
technical problems characterized by small Rm- 
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Figure 4: Illustration of the GIFT technique at a model of continuous casting, 
(a) Experimental set-up with the tundish, the single-port submerged entry nozzle, 
and the mould. The external field is produced by a coil, and the induced fields 
are measured at 7 fiuxgate sensors at each narrow face of the mould, (b) Four 
subsequent plots of the velocity field as reconstructed from the measured induced 
magnetic fields. The quantity A2, indicated by the grey scale, is the most appropriate 
measure for the identification of a vortex (Jeong and Hussain 1995). 
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